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Abstract 

We study the sensitivity of a long-base-line (LBL) experiment with neutrino 
beams from the High Intensity Proton Accelerator (HIPA), that delivers 10^^ POT 
per year, and a proposed IMt water-Cerenkov detector, Hyper-Kamiokande (HK) 
295km away from the HIPA, to the CP phase ((^^j^g) of the three-flavor lepton mix- 
ing matrix. We examine a combination of the narrow-band beam (NBB) at two 
different energies, {p-,^) = 2, 3GeV, and the NBB at {p-j^) = 2GeV. By allocating 
one year each for the two beams and four years for the 17^ beam, we can efficiently 
measure the — > z^g and transition probabilities, as well as the and 

17^ survival probabilities. CP violation in the lepton sector can be established at 
4(T (3it) level if the MSW large-mixing-angle scenario of the solar-neutrino deficit is 
realized, \5^J or \6^^^ - 180°| > 30°, and if 4\U,,\Hl - \UJ^) = sin^ 26^^^ > 0.03 
(0.01). The phase S-^^^ is more difficult to constrain by this experiment if there 
is little CP violation, 6^^^ ~ 0° or 180°, which can be distinguished at la level if 
sin2 20^^^ ^0.01. 
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Neutrino oscillation experiment is one of the most attractive experiments in the first 
quarter of 21st century. Many experiments will measure precisely the model parameters 
in the neutrino oscillations. In this article, we discuss the sensitivity of a long-base-line 
(LBL) experiment with conventional neutrino beams to measure the CP phase in the lepton 
sector. 

The Super-Kamiokande (SK) collaboration showed that the created in the atmo- 
sphere oscillates into z/^ with almost maximal mixing |^ . The SNO collaboration reported 
that the i/g's from the sun oscillate into the other active neutrinos [0. A consistent picture 
in the three active-neutrino framework is emerging. 

In the three-neutrino-model, neutrino oscillations depend on two mass-squared differ- 
ences, three mixing angles and one CP violating phase of the lepton-flavor mixing (Maki- 
Nakagawa-Sakata (MNS) |^) matrix. These parameters are constrained by the solar and 
atmospheric neutrino observations. One of the mixing angles and one of the mass-squared 
differences are constrained by the atmospheric-neutrino observation, which we may label 
0] as sin^ and 5m^^j^, respectively. The K2K experiment, the ongoing LBL neu- 
trino oscillation experiment from KEK to SK, constrains the same parameters Their 
findings are consistent with the maximal mixing, sin^ 29 ~ 1 (sin^ 9 ~ 0.5) and 

C) ATM ^ ATM ' 

'^"^ATM ~ (2 ~ 4) X 10~^(eV^). The solar-neutrino observations constrain another mix- 
ing angle and the other mass-squared difference, sin^ 26*^^^ and Sml^_^ , respectively. Four 

SOL SOL 

possible solutions to the solar-neutrino deficit problem are found: the MSW || large- 
mixing-angle (LMA) solution, the MSW small-mixing-angle (SMA) solution, the vacuum 
oscillation (VO) solution 0, and the MSW low-Sm'^ (LOW) solution. The SK collaboration 
0] and the SNO collaboration suggested that the MSW LMA solution is the most favor- 
able solution among them, for which sin^ 2^ = 0.7 ~ 0.9 and Smi = (3 ~ 15) x 10^^ 
eV^. For the third mixing angle, only the upper bound is obtained from the reactor neutrino 
experiments. CHOOZ [0 and Palo Verde found sin^ 2^^,^^ < 0.1 for 6ml^^ ~ 3 x lO'^ 
eV^. No constraint on the CP phase (5^j^g) has been reported. 

Several future LBL neutrino-oscillation experiments have been proposed to 

confirm the results of these experiments and to measure the neutrino oscillation parameters 
more precisely. One of those experiments proposed in Japan makes use of the beam from 
High Intensity Proton Accelerator (HIPA) |l^ and SK as the detector [^. The facility 
HIPA [0] has a 50 GeV proton accelerator to be completed by the year 2007 in the 
site of JAERI (Japan Atomic Energy Research Institute), as a joint project of KEK and 
JAERI. The proton beam of HIPA will deliver neutrino beams of sub-GeV to several GeV 
range, whose intensity will be two orders of magnitudes higher than that of the KEK 
PS beam for the K2K experiment. The HIPA-to-SK experiment with L=295 km base- 
line length and {E,y) ~ 1 GeV will measure at about 3 % accuracy and sin^ 9 at 
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about 1 % accuracy from the survival rate, while u -to-u^ oscillation can be discovered if 



sm^2e^^^ = 4|f/23|(l- |[/23|) ^0.006 pi. As a sequel to the HIPA-to-SK LBL experiment, 



prospects of using the HIPA beam for a very long base-line (VLBL) experiments with the 
base-line length of a few thousand km have been studied 0, ^7\, 0. Use of narrow-band 
high-energy neutrino beams {{E^) = 3 ~ 6GeV) and a lOOkton-level water Cerenkov 
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detector will allow us to distinguish the neutrino mass hierarchy (the sign of ml — 
if sin^ 2^^^^ ^0.03 Q]. If the LMA solution of the solar neutrino deficit is chosen by 



the nature, we can further constrain the allowed region of the S^^^ and sin 26'^^^^ |Q. 
However, because — > z7g appearance is strongly suppressed by the matter effect at such 
high energies, the measurement is not sensitive to the CP violating effects, ~ sin S 

o o ; o 7 MNS 

In this paper, we study the capability of an LBL experiment between HIPA and Hyper- 
Kamiokande (HK), a megaton- level water Cerenkov detector being proposed to be built at 
the Kamioka site |T^. Here a combination of the shorter distance (L = 295km) and low 
i/-energy {{Ej^) ~ IGeV) makes the matter effect small, and the comparison of z/^ — ^ z/^ 
and z7g appearance experiments is expected to have sensitivity to the CP violation 

effects proportional to sin S^^^ . 

The MNS matrix of the three-neutrino model is defined as 

3 3 
i=l i=l 

where a = e, fi,T are the lepton-flavor indices and u- {i = 1, 2, 3) denotes the neutrino mass- 
eigenstates. The 3x3 MNS matrix, Vmns; has three mixing angles and three phases in 
general for Majorana neutrinos. In the above parameterization, the two Majorana phases 
reside in the diagonal phase matrix V, and the matrix U, which has three mixing angles 
and one phase, can be parameterized in the same way as the CKM matrix [^. Because 



the present neutrino oscillation experiments constrain directly the elements, Ue2, f^es, and 
Ufj,3, we find it most convenient to adopt the parameterization |^ where these three 



matrix elements in the upper-right corner of the U matrix are chosen as the independent 
parameters. Without losing generality, we can take Ue2 and to be real and non-negative 
while f/e3 is a complex number. All the other matrix elements of the U are then determined 



by the unitary conditions 21 



The probability of finding the flavor-eigenstate f3 at base-line length L in the vacuum 
from the original flavor-eigenstate a is given by 



p 



where 
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satisfy A^^g + ^23 + ^31 = (^"^12 + ^"^23 + '^'^3i)(-^/2-E'jy) = 0. The two independent 
mass-squared differences are identified with the two "measured" ones, as follows; 



5m 
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(4) 



With the above identification, the MNS matrix elements are constrained by the observed 
survival probabilities, P^^^^y^ from the atmospheric neutrinos [Q, Py^^n^ from the reactor 
anti-neutrinos [llO| , and P^^^^^ from the solar neutrinos 0. The four independent 
parameters of the MNS matrix are then related to the observed oscillation amplitudes as 



\U, 



e3 



1 - sin^ 26^ 
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1 ± a/1 - sin^ 29, 
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arg {Ues 
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(l-|f^e3| 



sin^ 29. 
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(5a) 
(5b) 

(5c) 
(5d) 



The CP phase of the MNS matrix, S^^^, is not constrained. The solution eq.(^) follows 
from our convention [|], Uei > Ue2, which defines the mass-eigenstate u^. In this convention, 
there are four mass hierarchy cases corresponding to the sign of 6mfj] I {Sml^ > 5m\2 > 0), 
II {6ml^ > > dmj^), III {dmj^ > > dmj^), and IV (0 > dmj^ > dmj.^) §. If the MSW 
effect is relevant for the solar neutrino oscillation, then the neutrino mass hierarchy cases 
II and IV are not favored. When sin^ 26*^^^^^ 7^ 1, there is an additional twofold ambiguity 
in the determination of ?7^3 in eq.(|5BD. In order to avoid the ambiguity, we adopt the 
element itself, or equivalently sin^6'^^j^ defined in eq. (|5E|) , as an independent parameter 
of the MNS matrix. Summing up, we parametrize the three-flavor neutrino oscillation 
parameters in terms of the 5 observed (constrained) parameters Sm"^^^^, Sm"^^^, sin^^^^^^, 
sin^ 29 , sin^ 29 and one CP- violating phase 6 , for four hierarchy cases. 

sol' rot MNS' 

Neutrino-flavor oscillation inside of the matter is governed by the Schrodinger equation 
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where Hr, is the Hamiltonian in the vacuum and a is the matter effect term 



2V2Gpn^E^ = 7.56 x 10"^(eV2) 



g/cm^ 



\GeVj 



(7) 



Here is the electron density of the matter, E^^ is the neutrino energy, Gp is the Fermi 
constant, and p is the matter density. In our analysis, we assume for brevity that the 
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density of the earth's crust relevant for the LBL experiment, between HIPA and HK is a 
constant, p = 3, with an overall uncertainty of Ap = 0.1; 



p (g/cm^) = 3.0 ±0.1. 



The Hamiltonian is diagonalized as 
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by the MNS matrix in the matter U . The neutrino-flavor oscillation probabilities in the 
matter 



Pu 



'a3 



(10) 



takes the same form as those in the vacuum, with Ajj = (Aj — Xi)L/2E,j, if the matter 
density can be approximated by a constant throughout the base-line. Because the effective 
matter potential for anti-neutrinos has the opposite sign with the same magnitude, the 
total Hamiltonian H governing the anti-neutrino oscillation in the matter is obtained from 
H as follows [il , 



H (5ml^, 5ml^) = -H* {-6ml„ -6ml,) . (11) 

We make the following simple treatments in estimating the signals and the backgrounds 
in our analysis. 

• We assume a 1 Mega-ton water Cerenkov detector, which is capable of distinguishing 
between CC events and p^ CC events, but cannot distinguish their charges. 

• We do not require capability of the detector to reconstruct the neutrino energy. 

Although the water Cerenkov detector has the capability of measuring the energy of the 
produced p and e as well as a part of hadronic activities, we do not make use of these 
information in this analysis. We only use the total numbers of the produced p^ and 
events from z/^ or narrow-band-beams (NBB). The NBBs from HIPA deliver 10^^ 
protons on target (POT) in a typical 1 year operation, corresponding to about 100 days 
of operation with the design intensity |[TB|. Details of the NBB's used for this study are 



available from the web-page |23 



In the following discussion, we examine z/^ NBB's with the mean vr momentum (p^^) = 
2GeV (NBB(2GeV)) and {p^) = 3GeV (NBB(3GeV)), and NBB with (p^) = 2GeV, 
(NBB(2GeV)). For our input ('true') value of 5"^^^^ = 3.5 x 10"^ eV^, the probability 
Pu^^u, has a broad peak at ~ IGeV. NBB(2GeV) and NBB(2GeV) are chosen to 
maximize the transition probability, since (E^) ~ {pt,)/2. Because Pu does not change 
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NBB {{p^) = 2GeV) 


NBB {{p^) = 3GeV) 


NBB((p^) = 2GeV) 














0.06 


0° 


5.0 X 10^ 


8.5 X 10^ 


1.6 X 10^ 


1.1 X 10''^ 


1.6 X 10^ 


2.2 X 10^ 




90° 


5.1 X 10^ 


5.9 X 10^ 


1.6 X lO'' 


8.0 X 10^ 


1.6 X 10^ 


2.8 X 10^ 




180° 


5.1 X 10^ 


7.9 X 102 


1.6 X 10^ 


9.1 X 10^ 


1.6 X 10^ 


2.0 X 10^ 




270° 


5.1 X 10^ 


1.1 X 10^ 


1.6 X 10^ 


1.2 X 10^ 


1.6 X 10^ 


1.5 X 10^ 


0.01 


0° 


5.1 X 10^ 


1.7 X 10^ 


1.6 X lO'' 


2.3 X 10^ 


1.6 X 10^ 


4.5 X 10^ 




90° 


5.1 X 10^ 


6.2 X 10^ 


1.6 X 10^ 


9.5 X 10^ 


1.6 X 10^ 


6.7 X 10^ 




180° 


5.1 X 10^ 


1.4 X 10^ 


1.6 X 10^ 


1.4 X 10^ 


1.6 X 10^ 


3.7 X 10^ 




270° 


5.1 X 10^ 


2.5 X 10^ 


1.6 X 10^ 


2.7 X 10^ 


1.6 X 10^ 


1.6 X 10^ 



Table 1: Expected number of CC signal events from i/^ — > z/^, z^e oscillations for 
NBB(2GeV), NBB(3GeV) and those from 17^ ^ TT^.TT^ oscillations for NBB(2GeV), with 
IMt-year exposure. The results are shown for the parameters of eq.(0). 

much in the range E^^ ^ 0.6 ~ 1.2 GeV, our results do not depend strongly on the true 
value of the ^^^^tm- ^ong as it stays in the range (2 ~ 5) x 10~^ eV^ 

The signals in this analysis are the numbers of and CC events from NBB (2, 3GeV) 
and those of the and CC events from NBB(2GeV). These are calculated as 

/■lO GcV 

N,{v^- (pj) = MNa / dE,^,^{E,- {p^))P,^^,^{E,)a'^^{E,) , (12a) 



rlO GcV 

Nji^,, (Pn)) = MNa / dE^^^^{E^- {pJ)P^^^^^{E^)a^^{Err) , (12b) 



for / = e or /i, where M is the mass of detector (IMega-ton), Na = 6.017 x lO^^ is the 
Avogadro number, ^^^^{E^; (p^)) and $77^(^57; {Pn)) are the flux of in NBB((p^)GeV) 
and in NBB((p^)GeV), respectively. The flux is negligibly small at E^^ > lOGeV for 
the NBB's used in our analysis. The cross sections are obtained by assuming a pure water 
target g]. 

Typical numbers of expected CC signals are tabulated in Table |1] for the parameter 
setsQ: 

= ^Ktm = 3-5 X 10-3 eV' , (13a) 
= 5mJ^^ = 1.0 X 10-%V^ (13b) 
0°, 90°, 180°, 270°, (13c) 

p = 3g/cm^ (13d) 

* Recently KamLAND collaboration confirmed that only the LMA solution of the solar- neutrino deficit 
problem is consistent with the data||2^. The allowed region of Sm^^^ is found to be either (6 — 9) or 
(13 — 19) X lO^^eV^, slightly below or above our input value. The conclusions of this paper remain valid 
no matter which region its true value is. 



sin2 2^3^^ = 0.8, Sml,= 
sin^2e ^ = 0.06, 0.01, S , = 

RCT ' ' MNS 
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The numbers in the Table |l| are for 1 Mt-year exposure with 10^^ POT per year for 0.77 
MW operation of HIPA at L = 295 km. From Table |l|, we learn that the transition events, 
Ne and iVg, are sufficiently large to have the potential of distinguishing the CP conserved 
cases, 5 =0° and 180°, from the CP violating cases of 5 = 90° and 270°, even if 

' MNS ' ° MNS ' 

sin^26'j^^^ = 0.01. We also find that the survival events, A^^ and A^^, barely depend on 
the CP phase. The ratio iV^(2GeV)/A'^,(2GeV) is approximately a^^ /a^^ ~ 2.9, because 
both the fiux and the survival rates are approximately the same for and . From the 
comparison of Ni{2GeY) and Ni{2,GeY), we find that A'^(3GeV)/A'^(2GeV) ~ 3 because 
of the rise in the cross section (~ 1.5) and the increase in the survival rate (~ 2). The 
appearance signal A^^e increases only slightly at higher energies because a slight decrease 
in the transition probability cancels partially the effect of the rising cross section. Most 
notably, we find that the difference between the predictions of = 0° and 180° cases is 
significantly larger for Neiv^; {p^) = 3GeV) than that for Neiv^; {p^) = 2GeV). 

The above results can be seen clearly in Fig.|I|, where we show the expected number 
of z7e CC events Ng for NBB(2GeV) with 4Mton-year plotted against those of the CC 
event A^^ for NBB(2GeV) (left) and for NBB(3GeV) (right), both with IMton-year. The 
CP-phase dependence of the predictions are shown as closed circles for the parameters of 
eq.(|l3D at sin^ 2^^^^ = 0.06, 0.04, 0.02, and 0.01. Comparable numbers of 17^ CC events 
(A^e) and z/g CC events (A'e) are expected by giving 4 times more 17^ than z/^ beams. At 
each sin^ 26*^^^^ the events are expected to be smaller at ^^^^ = 90° (solid- squares) 

than at 6^^^ = 270° (open-squares). The trend is opposite for the z7^ — > z/g events, and 
thus anti-correlation allows us to distinguish the two cases clearly. On the other hand, 
the expected number at S^^^ = 0° (solid-circles) and that at 6^^^ = 180° (open-circles) 
do not differ much for NBB(2GeV) and NBB(2GeV). We find that NBB(3GeV) predicts 
significant differences between the two CP-invariant cases without loosing event numbers. 

In this report, we assume IMton-year exposure each with NBB(2GeV) and NBB(3GeV) 
and 4Mton-year exposure of NBB(2GeV), and examine the capability of HIPA-to-Hyper- 
Kamiokande experiments to measure the CP phase, S^^^, under the following simplified 
treatments of the backgrounds and systematic errors. 

For the z/^ and z/^ CC signal from NBB(z/^; (p^)), Ne{iy^; {pn)) and N^{iyf^; {pn)), respec- 
tively, we consider the following backgrounds: 

Ar,((pJ)BG = Ar,(z/,; (pj) + Ar^(z7^; (pj) + Ar^(z7^; (pj) + N^reiNC; (pj) , (14a) 
N,i{p^)ho = N,{v,- ip^)) + NjriV^- ip^)) + Ar^(z7,; (pj) . (14b) 

The first 3 terms in the r.h.s. are calculated as 

(\ rlO GeV 

iV(-)(^J; (p.)) = MNa / dE^^UEu. (P.))P(-) (-)(i?.)ap_^(i?,) , (15) 

I Jo Vc, 
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Figure 1: The CP phase dependence of Neiv^] 2GeV) for 4 Mt-year plotted against 
A^el^'/i! 2GeV) for 1 Mt-year in the left figure, and against iVe(t'^; 3GeV) for 1 Mt-year 
in the right figure. 5^^^ = 0° (solid-circle), 6^^^ = 90° (solid-square), 6^^^ = 180° (open- 
circle), and 6^^^ = 270° (open-square). The results are for the parameters at eq.(0). 

where and ^-p^ stands, respectively, for the secondary and flux of the primarily 
z/^ NBB. The last term in eq. ( |14a| ) for the e-like events gives the contribution of the NC 
events where produced vr^'s mimic the electron shower in the HK. By using the estimations 
from the K2K experiments we use 

with 

P,/^^ = 0.25 X (1± 0.1)%, (17) 

where the NC event numbers are calculated as in eg. ([T5|) by replacing a^j-^ by a^^. The 10% 
error in the misidentification probability of 0.25% is accounted for as a systematic error 
||r5| . The r-lepton contribution is found to be negligibly small for the NBB's considered in 
this analysis. The background for the enriched beam NBB(2GeV) are evaluated in the 
same way. 

Summing up, the event numbers for each energy neutrino and anti-neutrino NBB's are 
calculated from the sum : 

N^iip^)) = N,{v^- (pj) + iV,((pJ)BG , (18a) 
iV,((pJ) = iVy(77;; (pj) +iV,((pJ)BG- (18b) 

Most importantly, we do not require the capability of the HK detector to distinguish 
charges of electrons and muons. In Table |^ the expected numbers of CC and NC events 
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NBB((pJ) 












NBB(2GeV) 
IMton-year 


CC 


2.8x10^(1) 


2.2 X 10^(0.008) 


1.9 X 10^(0.007) 


1.3 X 10^(0.0005) 


NC 


1.1x10^(1) 


8.1 X 10^0.007) 


8.1 X 10^(0.007) 


5.3(0.0004) 


NBB(3GeV) 
IMton-year 


CC 


4.5x10^(1) 


3.1 X 102(0.007) 


2.0 X 10^(0.004) 


1.5 X 10^0.0003) 


NC 


1.6x10^(1) 


1.1 X 10^(0.006) 


8.6 X 10^(0.005) 


6.3(0.0004) 


NBB(2GeV) 
4Mton-year 


CC 


3.0 X 103(0.09) 


1.9 X 10^(0.005) 


3.5x10^(1) 


2.5 X 10^(0.007) 


NC 


1.2 X 10^(0.08) 


6.9 X 10^(0.005) 


1.5x10^(1) 


1.0 X 10^(0.007) 



Table 2: Expected number of the CC and NC events at HK in the absence of oscillations. 
The results are for 1 Mton-year for the enriched NBBs and 4Mton-year for enriched 
NBB from HIPA. The numbers in the parenthesis give the fraction of each mode against 
the main mode whose numbers are shown by bold letters. 

at HK in the absence of oscillations are shown for IMton-year each for the enriched 
NBB's and 4Mton-year for enriched NBB. The event numbers from the main (enriched) 
neutrinos are shown by the bold letters. The numbers in the parenthesis are the fractions 
as compared to the corresponding main mode. From the comparison between NBB(2GeV) 
and NBB(2GeV), we find that the fraction of the secondary-beam contributions is much 
larger for the z7^-beam than that for the z/^-beam. This is essentially because CC cross 
section is about a factor of three smaller than the CC cross section at Ej^ ~ IGeV. 

In Fig.|^, we show the expected — z/g {V^ Ve) signal and background event num- 
bers for the parameters of eq.(|l^) for sin^ 26'^^^^ = 0.01 ~ 0.06. The solid-circles show the 
number of expected signal events for 5^^^ = n x 10° (n = 1 ~ 36). The numbers of signal 
events are largest at around 6^^^ = 270° for NBB(2,3GeV), while those for NBB(2GeV) 
are largest at around 90°, as is expected from the CP phase dependence of A^^e and N-g 
shown in Fig.|l]. The open-triangle denotes CC events, which give the largest 

background for the experiments with NBB(2,3GeV), and the second largest background 
for NBB(2GeV). The open-square denotes z/g T^e CC events that gives the largest back- 
ground for NBB(2GeV), but is negligible for NBB(2,3GeV). The open-diamond denotes 
the background from the NC events, where vr^'s are miss-identified as electrons. They give 
the second largest background for NBB(2,3GeV). Backgrounds from —>■ Vf, transition 
events for NBB(2,3GeV) and those from z/^ v^. transition events for NBB(2GeV) are 
shown by open-circle. These transition backgrounds depend on the CP phase and they 
tend to cancel the ^^^^ dependence of the signals, but their magnitudes are small. The 
background level starts dominating the signal at sin^ 26*^^^ ^0.02. 

The background numbers for the /x-like signals are found to be negligibly small (~ 10^^) 
for NBB(2, 3GeV). Those for NBB(2GeV) are found to be about 21% of the signal almost 
independent of sin^ 2^ . In Both cases, the major background comes from the secondary 
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Figure 2: The sin^ "^^ncT dependence of the expected signal and background event numbers 



for the parameters of eq. (|T3|) for sin 26 
number of expected signal events for 6 



MNS 



= 0.01 ~ 0.06. Solid-circles stand for the 
n X 10° (n = 1 ~ 36). Open-diamonds denote 



the 71 background from the NC events. Open-triangles and open-square show Ue and z7g 
survival events. Open-circles are z7^ Ve transition events for NBB(2, 3GeV) and — > z/g 
transition events for NBB(2GeV). 



i^fi) survival events. 

Our analysis proceeds as follows. For a given set of the model parameters, we calculate 
the expected numbers of all the signal and background events for each NBB((p7r)) and 
NBB((p7r)), by assuming 100% detection efficiencies for simplicity. The resulting numbers 
of yU-like and e-like events are then denoted by Nj["^^{{p.„)) and -/V'g™'^((p^)) for NBB((p^)), 
and iV;:*™'((p.)) and K'"'"' {{p.)) for NBB((p.)). 

We account for the following two effects as major parts of the systematic uncertainty 
in this analysis. One is the uncertainty in the total flux of each neutrino beam, for which 
we assign the uncertainty, 

/l„((P.)) = l±0.03, (19) 

independently for = u^, iy^,V^,V^ and for NBB(2GeV), NBB(3GeV), and NBB(2GeV). 
Although it is likely that correlation exists among the flux uncertainties, we ignore pos- 
sible effects of correlations in this analysis. By using the above flux factors, theoretical 
predictions for the event numbers, N[^^{{p^)) and A^f'*((p^)), are calculated as 

Ntiip^)) = Ui{pJ)N,{u,, (p.)) + fuS{p^))Ne{u^, (p^)) 

+f^^{{p^))Njiv,, (pj) + f^^i{p,))Njiv^, (pj) 

+ke Pe/NC E /.J(p.))iV.7((p.)) , (20a) 
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+5.,e Pe/NC E7..((P.))<^((P.)) ' (20b) 

where the last terms proportional to 6i^e are counted only for i = e. As the second 
major systematic error, we allocate 3.3% overall uncertainty in the matter density along 
the base-line, eq.(P). The fit functions are hence calculated for an arbitrary set of the 6 
model parameters, the 12 flux normalization factors, and the matter density p. 
The function of the fit in this analysis can now be expressed as 




2 

'fuA{Pn)) 



+ E 



NBB 




+ E 

Ki(p ^))-K'i(p^)) y I ( Ni'\{p.))-Nr i{p.)) 

7. ((p.)) - 1.0^ 



E 



0.03 



/p_3.0V /5m2 -5m2 /sin2 2e-^^* -sin2 2^*™n^ 

/ y tj.w I / SOL SOL | i | SOL SOL | (21) 

V 0.1 / V 0.1 X Wi™" / I 0.06 / ' ^ ' 



SOL 



where the summation is over NBB(2GeV), NBB(3GeV) and NBB(2GeV). Even though 



we have only one NBB in our analysis, we retain the summation symbol in eg. (pTD for 
the sake of clarity. The last two terms are added because KamLAND experiment ||2^ will 



measure Sm'l at 10% level and the solar neutrino experiments constrain sin^ 26' with 



the la error of about 0.06 for the LMA parameters of eg. (p!3D . The individual error for 
each N^{{p^)) {N^{{p^))) is statistical only, whereas the error for each Ne{{p^)) (A^e((P7r))) 
is a sum of the statistical errors and the systematic error coming from the 10% uncertainty 
in the e/vr^ misidentification probability of eg. (p!7|) , 



a,{{p^)) = ^Nf[-^{{p^)), (22a) 
ae ( (p.) ) = ^Arf-((K)) + {o.lNt^riNC; {p^))f . (22b) 



The errors for the NBB(2GeV) case are calculated similarly as above. 

We show in Fig.|^ regions allowed by the HIPA-to-HK experiment in the plain of 
sin^ 29 ^ and 6 ^ . The mean values of the input data are calculated for the LMA pa- 

RCT MNS ^ 

rameters of eg. ([131). In each figure, the input parameter point (sin^ 2^*™^ ''^mnr) shown 
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Figure 3: Regions allowed by the HIPA-to-HK experiment are shown in the plain 
of sin^ and ^j^^^g- The assumed experimental conditions are 1 Mt-year each for 

NBB(2GeV) and NBB(3GeV), and 4 Mt-year for NBB(2GeV) with lO^i POT/year. The 
input data are calculated for the LMA parameters of eq.(|TB|). In each figure, the input 
parameter point (sin^ 26'^™^ , 5^^g) is shown by a solid-circle for sin^ 2^^™^ = 0.06, and by 
a solid-square for sin^ 26'^™^ = 0.01. The regions where Xmin 4, and 9 are depicted by 
solid, dashed, and dotted boundaries, respectively. 
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Figure 4: The same as Fig.|^ but with 2 Mt-year for NBB(2GeV) and 4 Mt-year for 
NBB(3GeV) only. The two- fold ambiguity in the fit is clearly seen. 

by a solid-circle for sin^ 2^*™^ = 0.06, and by a solid-square for sin^2^*™*^ = 0.01. The 

RCT ' ^ RCT 

regions where xLin 4, and 9 are depicted by solid, dashed, and dotted boundaries, 

respectively. All the 6 parameters, 5"^^™*' ^^^^^atm' '^"^sol**' ^^^^^^sol ^^^^^^rct' ^mns' 
the matter density p-^**, and the 12 flux normalization factors are allowed to vary freely in 
the fit. 

From the top-right and bottom-right figures for 5^^^ = 90° and 270° respectively, we 
learn that S^^^ can be constrained to ±30° (±60°) at the la (Scr) level, even if sin^ 26'^™^ = 
0.01. This is because A^e ± Ne constrain sin^ 26*^^^^ and Ne/Ne distinguishes between 
6 = 90° and 270° in Fig.|3, whereas the remaining parameters (Sm'^ and sin^ 6 ) are 

MNS t> a; t> Jr' \ ATM ATM ' 

constrained by the and survival data, A^"^ and A^^. The accuracy of the S^^^ mea- 
surement does not decrease significantly for sin^2^^™^ = 0.01 despite the large background 
level, because the 5^^^^^ -dependence of the signal exceeds significantly the 3% uncertainty of 
the background level from the fiux normalization factors in eg. ([T9|) . We find that the CP 
violation signal can be distinguished from the CP-conserving cases {S^^^ = 0° or 180°) at 
4a (3a) level for all 5 values in the region \5 1,15 - 180° I > 30° if sin^ 2^*'""'^ ^0.03 

^ ' MNS ° I MNS I ' I MNS ' RCT ~ 

(0.01), for the LMA parameters of eq.(0) and for the systematic errors assumed in this 
analysis. 

The situation is quite different for the CP-conserving cases of 5*'^"'^ = 0° or 180° shown 

^ ° MNS 

in the left-hand side of Fig.|[ ^^j^g can be constrained to better than ±7° (11°) accuracy 
at la level for sin^ 26*^^^^ ^0.06 (0.01), but the two cases cannot be distinguished at 2a 
level. This is mainly because of the similarity of N^/Ne between S^^^ = 0° and 180° in 
Fig.|I|. The difference between the two cases is larger for NBB(3GeV). If we remove the 
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NBB(3GeV) data from the fit, we find that the two cases cannot be distinguished even at 
la level. This two-fold ambiguity between 6 and 180° — 6 is found in general for 

MNS MNS 

all 5j^j^g, because the difference in the predictions can be adjusted by a shift in the fitted 
sin^ 29^^^ value; see Fig.|l|. 

As a demonstration of the effect of using two NBB's, NBB(2GeV) and NBB(3GeV), 
in the analysis, we show in Fig.^ the fit results when the data are generated by using 
NBB(2GeV) and NBB(2GeV) only, each at 2 Mton-year and 4 Mton-year, respectively. It 
is clearly seen from the figures that the 'mirror' solution at S^^^ = 180° (0°) can fit the 
data as well as the 'true' solution at 6^^^ = 0° (180°). Essentially the same results are 
obtained when we replace NBB(2GeV) by NBB(3GeV) in the above analysis. It is only by 
combining the two NBB's that we can distinguish the two solutions as shown in Fig.^ We 
find less significant difference from the results of Fig.|3| when the input S value is 90° or 

o o l_l r- MNS 

270°. 

It is remarkable that the la error of 6 is as large as 30° for 5*™*^ = 90° and 270° while 

MNS ^ MNS 

it is less than 10° for 5*'""^ = 0° and 180°. This is simply because the S dependence of 

MNS •' MNS ^ 

the Ufj^-to-Ue (and also z7^-to-z7e) oscillation probability is roughly proportional to sin^^^j^^, 
in the vicinity of the first dip of the z/^-to-i/^ survival probability. 

We close this article by pointing out that the low-energy LBL experiment like HIPA- 
to-HK cannot distinguish between the neutrino-mass hierarchy cases (between I and III) 
because of the small matter effect at low energies. If we repeat the analysis by using the 
same input data but assuming the hierarchy III in the analysis, we obtain another excellent 
fit to all the data where the fitted model parameters are slightly shifted from their true 
(input) values. VLBL experiments at higher energies at L > 1000km @] are needed to 
determine the mass hierarchy. 
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